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Abstract 

A general result is presented on the lack of second order corrections and on the form of the leading 
order Floquet quasi-energies in the high-frequency approximation for a two-level driven system. 
Then, by a perturbative approach in the high-frequency approximation that uses dual Dyson series 
and renormalization group techniques [M. Frasca, Phys. Rev. B 68, 165315 (2003)], we obtain 
a third order correction for a sinusoidal driving field that is the same obtained by Barata and 
Wreszinski [Phys. Rev. Lett. 84, 2112 (2000)], confirming their result that localization deviates 
from the zeros of the zero-th order Bessel function when higher order terms are taken into account. 
Three different expressions for this correction are compared. An important consequence of this 
result is that gives complete support to the correctness of both methods. Finally, the limitation of 
high-frequency approximations is presented by comparison with numerical results for a sinusoidal 
driving field. 

PACS numbers: 72.20.Ht, 73.40.Gk, 03.67.Lx 
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In a recent paper by Barata and Wreszinski (BW)i a new method to solve perturbatively 
the driven two-level system was presented. These authors showed that a third order correc- 
tion, for the case of a sinusoidal driving field, appears to change the leading order Floquet 
quasi energies, in the limit of large frequency, in such a way to modify the naive result that 
one has localization in the two-level system at the zeros of the zero-th order Bessel function. 

This analytical result is in agreement, to a certain extent as we will see, with a work by 
Crefneld about the crossing manifolds of the localization points in this model 2,3 . A different 
work by Platero, Brandes and Aguado^ showed that analytical and numerical results agree 
fairly well taking into account the result by Barata and Wreszinski. 

The question of higher order corrections is a relevant one in view of the needs of workable 
two-level systems to be manipulated as qubit in quantum computers. This means that 
general results are important to know the proper way to drive a two-level system for such 
aims. 

In this paper we present a several results on higher order corrections to Floquet quasi- 
energies for a driven two-level system. The approach we use has been already presented irA 
The method given there uses the dual Dyson series and renormalization group techniques. 
We improve the computations till third order showing equivalence with the method of Barata 
and Wreszinski giving in this way fully support to both methods in the high-frequency 
approximation. 

A two-level driven system that is generally considered to implement a qubit by a solid 
state device has the form: 



being A the separation between the two levels, g the coupling constant and f(t) = f(t + T) 
a periodical driving field with period T and to = 2n/T and ai,a^ Pauli matrices. The high- 
frequency regime is identified by the condition - < 1. Floquet theory can be applied and 
from this we know that localization appears at the crossing of the quasi-energies^. Then, a 
proper choice of the parameters can permit to drive a qubit. The other crucial parameter 
for our aims is z — 2g/u that we will use in the following. 

We can prove a general result on the model (P) using the method of dual Dyson series 
given in Ref.— . We introduce the unitary transformation 



H — —cr 3 + ga\f(t) 
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that gives the transformed Hamiltonian 



A 

H F = —a 3 exp 



-2giax f dt'f{t' 
Jo 



(3) 



By introducing the Fourier series 
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where U\ can be treated as a c-number, it is not difficult to write down the Hamiltonian (jSJ) 
in the form 

H F = H F0 + <x 3 £ a n (z)e- m ^ +a 2 ^ b n {z)e~ m ^\ (6) 

being H F0 = a 3 a (z) + a 2 b (z) = cr 3 c (z) the leading order Hamiltonian whose eigenvalues 
give the leading order Floquet quasi-energies. This expression is obtained by considering 
that the coefficients of the Fouries series depend on o\. Then, we are able to derive an 
estimation of the Floquet quasi-energies through eq.( |5| ) for any kind of driven two-level 
model. 

We can extend the above result by computing the second-order correction to the quasi- 
energies by the method of dual Dyson series^. The computation with Hamiltonian (jUJ) is 
quite straightforward giving the following secular contribution to the unitary evolution 



U s F ec (t,t ) = -i[a 3 a (z) + a 2 b (z)](t - t ) + ia x £ "M^M (, /n) 



nuo 
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where we can recognize that a second order correction cannot appear when a n = ±a_ n and 
b n = ±6_„. This is precisely the case for a sinusoidal driving field having c n (z) = J n (z) 
being J n (z) the Bessel functions of integer order. This means that eq.(29) and eq.(30) in 
Ref.— have not the second order correction and should be corrected by eliminating it. The 
same is true for a square wave driving field^ f(t) = O(t) — 2@(t — T/2) in the interval 
< t < T. In this case we have the Fourier coefficients 

sin7T2; ( — l) n COS7T2; — 1 



z 
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z 2 — n 2 



+ ^c^l- 



z 2 — n 2 
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and the leading order Hamiltonian H F0 = e (a 3 cos(ttz/2) + <j\ sm(irz/2)), being e 
— sin(7rz/2) and ±e the leading order Floquet quasi-energies. 



T 3 =i^r<Tz fdh f 1 dt 2 f t2 dt 3 ^^Y,M^WU^~ iUaximtl ~ nt2+P ^ (11) 



A sinusoidal driving field has f(t) = cos(ut). Floquet theory^ gives localization, as a 
first approximation, for Jq(z) = 0. One may ask if this result is true at any order or just as 
a first order approximation. Barata and Wreszinski obtained that a third order correction 
given by 

BW( ) = V ^ 2rt i+ 1 ( Z )^ 2 "2+l( 2: )^-2(ni+n2+l)(^) /q\ 

{) ~ ^\i% 2 (2ni + l)(2na + l) lJ 

that is not generally zero at the zeros of Jo (z). In a successive work Barata and Cortez^ 
proved that the same correction could be written as 

r>r<(~\ ^ 3 Jni( z )Jni-ri2( z )Jn 2 ( z ) / ln \ 

W = "4^ J%o ^ ' ( } 

Our method given in Ref.— needs the computation of the third order term of the unitary 
evolution by dual Dyson series 

/ dh / dt 2 / 

/to Jto Jto m n p 

starting from a generic initial time to and being J n Bessel functions of integer order. Without 
repeating all the procedure here, our approach gives the third order correction 

pf \ = A 3 J (z) Jn( z ) A 3 J ni (z) J ni +n 2 ( z ) Jn 2 ( z ) 

4cj 2 ^ n 2 8oo 2 , n ^! , ni n 2 ' 1 ' 

that yields the complete term with respect to eq.(31) in Ref.—. We have given explicitly 
the first term with Jq(z) that recovers the small z limit. This limit is recovered by all the 
expressions BW(z), BC(z) and F(z) as it should be. As already pointed out by Barata and 
Corteai, it is not possible to prove directly equivalence of this formulas by analytical means. 
The relationship does not depend on the properties of the Bessel functions. So, we need to 
content ourselves of a direct numerical check here. 

In fig^we have compared the three expressions above, normalized by the factor A 
and the result is absolutely satisfactory. 

So, the BW method and our method give both the same result supporting each other 
with respect to correctness as should be expected. To use one or the other becomes just 
a matter of taste. The fact that two different analytical approaches, in the limit of large 
frequency, give the same result permits us to conclude that, also by this way, it is definitely 
proved that zeros of Jq(z) are just approximately localization points . 
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FIG. 1: (Color online) The three terms, normalized by the factor — 1 are given here in red 
(solid) for F(z), green (points) for BC(z) and blue (dotted) for BW(z). The red and green curves 
are almost identical. 

One may ask how far this third order corrections can be pushed in the ratio A/a;. More 
generally, it is rather interesting to understand the limitations of such computations. The 
answer can be obtained from fig.© where is seen that for values of the ratio A/a; w 0.5 all 
the high-frequency approximations begin to fail demanding for higher order terms. 

In order to give a proper picture of the problems encountered in Ref.— we note as the series 
we obtained in it for the Floquet quasi-energies till third order was not correct [eqs.(29),(30) 
and (31) in Ref.—]. Rather we have shown here that no second order correction appears and 
that the third order correction deviates from the naive leading order result of localization 
at the zeros of the zero-th order Bessel function. These corrections make our approach 
completely equivalent, for the results, to the method conceived by Barata and Wreszinski 1 . 
It should be said that our method makes the computation of Floquet eigenstates and quasi- 
energies a completely algorithmic matter, relying all the calculations on a very simple recipe. 

In conclusion we have given general results for a driven two-level model giving a method to 
compute high-frequency leading order Floquet quasi-energies. Besides, we have proved fully 
equivalence between Barata and Wreszinski and our method to compute a high-frequency 
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FIG. 2: (Color online) The three terms, F in red (dotted), BC in green (solid) and BW in blue 
(diamonds) are compared with the numerical results on crossing manifolds in black (circles) given 
in Ref£. 



perturbation series. Finally, the limitation of these methods has been pointed out by com- 
paring them with a numerical computation. 
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